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Abstract 

We study the impact of next-to-next-to-leading order (NNLO) QCD corrections on par- 
tial decay rates in i? ^ X^^vi decays, at leading-order in the l/mf, expansion for shape- 
function kinematics. These corrections are implemented within a modified form of the 
BLNP framework, which allows for arbitrary variations of the jet scale /Xj ~ 1.5 GeV. Our 
analysis includes a detailed comparison between resummed and fixed-order perturbation 
theory, and between the complete NNLO results and those obtained in the large-/3o ap- 
proximation. For the default choice /ij = 1.5 GeV used in current extractions of \Vub\ 
within the BLNP framework, the NNLO corrections induce significant downward shifts 
in the central values of partial decay rates with cuts on the hadronic variable P+, the 
hadronic invariant mass, and the lepton energy. At the same time, perturbative uncer- 
tainties are reduced, especially those at the jet scale, which are the dominant ones at 
next-to-leading order (NLO). For higher values of /Xj and in fixed-order perturbation the- 
ory, the shifts between NLO and NNLO are more moderate. We combine our new results 
with known power-suppressed terms in order to illustrate the implications of our analysis 
on the determination of \Vub\ from inclusive decays. 



1 Introduction 



The CKM element \Vub\ is a fundamental parameter of flavor physics. It measures the 
strength of 6 — > m quark transitions and determines the length of the side of the uni- 
tarity triangle opposite to the angle (3. The combined work of many theorists and ex- 
perimentalists has allowed to determine this parameter from data on both inclusive and 
exclusive semi-leptonic b ^ u decays; a summary of current results can be found in, for 
instance, [1-3]. A characteristic feature of these analyses is that the value of \ Vub\ deduced 
from inclusive decays is consistently higher than that from exclusive decays. The need to 
resolve this discrepancy motivates systematic improvements on all fronts. 

The goal of this paper is to improve the theory predictions for the inclusive decays B 
X^ii^e by including recently calculated next-to-next-to-leading order (NNLO) perturbative 
corrections to the partial decay rates used in the extraction of \Vub\- In general, these 
partial rates are available only in the portion of phase space referred to as the shape- 
function region, where the hadronic final state is a jet carrying an energy of order rrih 
and an invariant mass squared of order mbAqcD- A local operator product expansion is 
not valid in this region, and results are obtained through a factorization formalism, as a 
convolution of perturbative hard-scattering kernels with non-perturbative shape functions 
[4-6]. Different approaches to this factorization have been put forth in the literature, going 
under the names of BLNP [7,8], GGOU [9], and the dressed-gluon exponentiation [10]. 

In this paper we implement the NNLO perturbative corrections within the BLNP 
framework. In this approach, techniques from soft-collinear effective theory (SCET) [11- 
13] are used to obtain an arbitrary partial decay rate as an expansion in l/m^, according 
to 

r. = rw + [rw + ri^) + ...]. (i) 

The first term is of leading order in the heavy-quark expansion, whereas the terms in 
the square brackets account for power-suppressed effects and are estimated up to order 
1/m^. The NNLO corrections studied in the current paper affect only the leading-order 
term Tu^ . For kinematic cuts limited to the shape-function region, this term obeys a 
factorization formula of the schematic form [7, 14-16] 

^^^^-[H-J]M®SM, (2) 

where the symbol (8> denotes a convolution. The factorization formula contains a hard 
function H, related to physics at the hard scale fih ~ "^t, a jet function J, related to 
physics at the intermediate scale /Xj ~ (mbAqcD)^^^ ~ 1-5 GeV, and a non-perturbative 
shape function S, describing the internal soft dynamics of the B meson. A thorough 
analysis of this leading-power term at next-to-leading order (NLO) in renormalization- 
group (RG) improved perturbation theory was performed in [8]. We can extend this 
analysis to NNLO by putting together a number of pieces, which we describe in Section 2. 
In addition to calculating the higher-order perturbative corrections, we also modify the 
BLNP framework to allow for variations of the arbitrary matching scale /ij at which the jet 
function is calculated. This allows us to study the perturbative uncertainties associated 
with this scale, and makes for a straightforward matching with fixed-order perturbation 
theory, where logarithms between the hard and intermediate scales are not resummed. In 
the numerical analysis in Section 3 we apply our new results to partial decay rates with 
cuts on the hadronic variable P+, the hadronic invariant mass, and the lepton energy. 
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Our main findings are that the dependence on the scale /ij is sizeable at NLO and still 
significant even at NNLO, and that for the default choice /ij = 1.5 GeV typically used in 
the BNLP framework the NNLO corrections tend to shift the partial decay rates downward 
by a significant amount. We also study some qualitative aspects of the perturbative series 
by comparing results obtained in resummed and fixed-order perturbation theory, and 
with those obtained in the large-/3o approximation. To illustrate the significance of our 
results for the extraction of \Vub\, in Section 4 we combine the NNLO corrections to 
the leading-order term with known power corrections and experimental data in order 
to extract sample values of \Vub\ for several partial rates. For the particular choice of 
intermediate scale /i, = 2.0 GeV, the effect of the NNLO corrections is to raise the central 
value of \Vub\ by slightly less than 10% compared to the results at NLO; for higher choices 
of the intermediate scale and in fixed-order perturbation theory, the corrections are more 
moderate. We summarize our findings in Section 5. 



2 Partial decay rates in SCET 

In this section we briefiy review the BLNP formalism as applied to inclusive semi-leptonic 
h ^ u decays. We begin by recalling the master formula for an arbitrary partial decay 
rate derived in [8], expressed in terms of the hadronic variables 



It is given by 



P_ = Ex + \Px\, P+ = Ex-\Px\- (3) 

max 



dP+ \T^{y^)+T^- y^ax>l/o. 



(4) 



where 



(2/.) = {Ms - P^f I dyy' [(3 - 2y) h + 6(1 - y) h + V h] 



2 ry 







n2 IT/ 12 /-ymax 

rf = (M^ - dyyo (5) 



X [(6y(l + yo) - Qy^ - Z/o(3 + 2yo)) /i + 6y(l -y)f2 + Z/o(3y - 2yo) /s] • 
The dimensionless variables y, yo, and ?/max are defined as 

_ p_ - p+ _ p"'^^ -p+ _ pr"" - ^+ _ 1 

y~ Mb-P+' Mb-P+ ' " Mb-P+ ~ Mb-P+' ^ ' 

and we have introduced the quantity Eq, which is the minimum allowed lepton energy. 
The values of |/max and the integration range for P+ depend on the specifics of the partial 
decay rate under consideration. For cuts P+ < Ap and on the lepton energy, we have 
Z/max = 1 and < P+ < min(Ap, Mb — 2Eq). For a cut on the hadronic invariant mass, 
Mx < Mq, we have 

min(M^,MoVP+)-P_, 
1/max- Mb-P+ 

and < P+ < Mq. 
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The scalar functions fi are obtained as an expansion in l/rrih. The leading-order term 
has the form (2) and reads 



fr^\P+,y) = Huiiy,mb,fif) dujymbJ{ymb{P+-uj),fif)S{u;,fif). (8) 

As written, (8) achieves a factorization of perturbative and non-perturbative physics. The 
hard functions H^i contain physics at the scale m^, the jet function J contains physics 
at the jet scale (mhAqcD)^''^, and the shape- function 5 is a non-perturbative matrix 
element in heavy-quark effective theory (HQET) [4,6]. In SCET, one usually assumes 
the parametric limit where the hard scale is much larger than the jet scale, in which case 
any choice of /i/ leads to large perturbative logarithms in either H or J. To solve this 
problem one first calculates these functions at matching scales where they contain no 
large logarithms, and then evolves them to a common scale /i/ using the renormalization 
group (RG). Also, within the effective-theory framework, it is natural to extract the non- 
perturbative shape function at a low scale /zq, and then evolve it to the scale fif. The 
RG evolution can be achieved using results from [7, 17] for the hard and shape functions, 
and [18] for the jet function. For convenience, we list the solutions here: 



^-2ar(M.,M/)exp 



2S{fih,fJ'f) - 2ar{fih,fJ'f) In 2ay(/i/„/i/) 



xHui{y,mb,fih) 



(9) 



J{p'^,f^f) = exp [-4^(/ii,/x/) + 2ay(/ii,/i/)] j(<9^^,/ii 



1 f p 



2 \ VJl 



-IeVj 



. nvj) 



5(cj,/i/) = exp [25(/io,/i/) + 2ay_y(yUo,Ai/)l — — ^ / du 



s) Jo 



The RG exponents read [7] 



^^rcusp(«) 



/3(a) J (3{a') 



da' 



da 



Tcusp (Q^) 



(10) 



and similarly for ay (a^j), but with Fcusp replaced by the anomalous dimensions of the 
hard (jet) function. We have defined rjj = 2ar{fii, fif), rjs = 2ar(/i/, /io), and j is the 
Laplace transform of J. Finally, the star distribution is defined as 



dp^ 



1 f p 



p^ \IJ? 



J2 \ V 



f{p') = I dp 

'0 



2/(p')-/(o) fp'y , /(o) fQ'^' 



Inserting the above results into the factorization formula (8), and using the general rela- 
tions [19] 



S'(/ii,/i2) + S'(/i2,/i3) = S'(/ii,/i3) -Fin — ar(/i2,/U3) 



'121 
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one is left with 

25 (/ih, yUj) - 25" (/ii, /io) + 2ay (/ii, /iq) - 2ay (/i^, /xo) 

xi7„i(y,m,,/i,)z/-"^^^'^'''^'^ (13) 

xj In + 9rj, f^i ] ^^^^ duj — -[ S{uj,fio), 

\ Hi 'J T{r]) Jo [P+-UJ \ fii J \^ 

where now t] = 2ar(/ii, /io)- Note that the dependence on the arbitrary factorization scale 
Hf has disappeared. 

Equations (4) and (13) are the master formulas for calculating the leading-power 
contribution to a given partial decay rate in RG-improved perturbation theory. The results 
are formally independent of the matching scales fih and /ij, but a residual dependence 
remains when truncating the perturbative expansion at a given order. Moreover, by 
setting Hh = Hi = Hi one recovers fixed-order perturbation theory directly from the 
resummed results. The product of matching functions [H ■ J) (fi) is combined into a single 
coefficient function C{fi), and the RG exponents serve to evolve the shape- function from 
the scale fio to /i. We shall discuss and compare the perturbative uncertainties in both 
resummed and fixed-order perturbation theory at LO, NLO, and NNLO in the following 
section. The scale ho plays a special role in the analysis. We shall choose this as the scale 
at which we model the non-perturbative shape function; more details will be given below. 

We can evaluate the master formula (13) at NNLO in RG-improved perturbation 
theory by gathering together a number of results available in the literature. By NNLO, 
we mean the approximation which captures all of the order terms in both the matching 
functions and the RG exponents. For the matching functions H^i and j this counting is 
unambiguous: they are both needed at two loops. The H^i to this order can be derived 
from the calculations in [20-23], and j was calculated to NNLO in [24]; for the convenience 
of the reader, we list the results in the Appendix. As for the RG exponents, to define 
their NNLO expansion requires assumptions about the matching scales Hh, f^i, and ho- 
By default, we shall assume the hierarchy Hh ^ l^i ^ f^o- In that case, to account for 
all of the pieces in the RG exponents requires r^usp to four loops and the anomalous 
dimensions a^j and ay to three loops. However, both 7' and Fcusp are known to one loop 
lower, which adds a small uncertainty to the analysis. For the missing pieces, we shall use 
the [1,1] Fade approximation described in the Appendix. In practice, we shall always use 
Ho = 1-5 GeV in next section, so that Hi ~ A^o- In that case the general expression (13) 
resums some higher-order logarithms which are not large. 

To evaluate the partial decay rates requires a model for the shape function S{lj,Ho)- 
Experimental information on the shape function is provided by data on the photon energy 
spectrum in B ^ Xg'j decays, and can be used to guide the functional form of the model. 
In addition, model-independent constraints are provided by the fact that moments of the 
shape function, defined as [7, 8] 

PUJO 

Mn{luo,h)= / duJUJ^S{uj,H)^ (14) 
Jo 

can be calculated in a local heavy-quark expansion, as long as ljo ^ Aqcd [7,25]. In the 
renormalization scheme referred to as the shape-function scheme [7,25], the shape- function 



-2ar(/ih,/ii)ln — 
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moments are used to define the heavy-quark parameters order-by-order in perturbation 
theory. Limiting ourselves to the first two moments, which is sufficient to order in 
the heavy-quark expansion, one has 



Mi(/i/ + A(/i/,/x),/x) 
Mo(^/ + A(/i/,/i),/i) 

M2(/i/ + A(/i/,At),At) 



A(/i/,/i) , 



(15) 



where A = MB—rrib and fi^ is related to the kinetic-energy parameter in HQET. Requiring 
that a given shape-function model correctly reproduces the moment relations (15) puts 
constraints on its parameters. These constraints depend on the order in perturbation 
theory at which the moments are evaluated. In the numerical analysis, we shall always 
use the constraints obtained from the two-loop moments. The moments to this order 
can be calculated as explained in [25], using the two- loop expression for the quantity 
s{L,fi) obtained in [26]. In addition to the perturbative expressions, one needs numerical 
values for A{fif, fi) and fi^{fif, fi). These can be determined from global fits for the HQET 
parameters in other renormalization schemes using perturbative conversion relations. For 
the choice of scales /ij = /i = /i*, the connection with the pole scheme reads [25] 



m 



pole 



1 + 



271 7tt^ 



36 36 



17 , 47 
12^^ 54 ^ 



+ 



A** 



7r 



13 
81 



17 , 10 
27^ 243 ^ 



2as(/i*) , / «s(/U^) 



19 

18 ~ ^ 



17 



9 



TT 



20 



17 , 14 
6 ^ 27 ^ 



In the rest of the paper, we will use the notation mi,{ii^,^^) 



ml 



and fil{fi^,fi^ 



(16) 



for the HQET parameters in the shape- function scheme, evaluated at the scale /i* = 
1.5 GeV. Recent HFAG numbers for these parameters are ml = (4.707^0 053) GeV and 
/i*^ = (0.216lo o7g) GeV^ [1]. They are determined from information on S — > XJui 
moments alone, under the assumption that moments from B Xg'-f decays should not 
be used in the global fits, on grounds that the measurements are typically made at values 
of the photon energy where shape-function effects are expected to be non- negligible [3]. 
To obtain numerical values for the non-diagonal parameters mb{fif, /i) and fi^fif, one 
uses the generalization of (16) given in [25], along with the fact that the pole scheme 
parameters are scale independent. 

In the analysis that follows, we model the shape function as 



exp — 7- 



buj 
A 



(17) 



The normalization factor A/", as well as the model parameters b and A, are tuned to satisfy 
the two-loop moment constraints described above. To get a feeling for how this compares 
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Figure 1: Left: model of the shape-function with the moment relations evaluated at one- 
loop order (dashed) and two-loop order (solid), for ml = 4.71 GeV. Right: model of the 
shape-function with the moment relations evaluated at two loops, for ml = 4.66 GeV 
(dotted), ml = 4.71 GeV (sohd), and ml = 4.77 GeV (dashed). 



with the model generated from the one-loop moments and used in [8], we show in the 
left-hand plot of Figure 1 the shape- function model (17) tuned to reproduce the moment 
relations for ml = 4.71 GeV and /x*^ = 0.2 GeV^ at one- and two-loop order. In both 
cases, we use the four-loop running coupling with ^^(M^) = 0.1176, and match onto the 
four-flavor theory at 4.25 GeV. The results also depend on the UV cutoff of the moment 
integration range, which we choose as uo = Mb—2Eo, with Eq = 1.8 GeV. We furthermore 
take /io = 1.5 GeV. The figure shows that the shape- function model generated using the 
two- loop moment relations has a smaller average value for u) < 1.0 GeV compared to that 
generated with the one-loop moments. Since to leading order in a given partial rate is 
directly proportional to the integral of the shape function over a window specified by the 
cut, the model with the two- loop constraints tends to yield lower values for partial decay 
rates with kinematics restricted to the shape-function region. In the right-hand plot of 
Figure 1 we show the shape-function model generated using the NNLO moment relations 
for three different values of ml, which cover the range of the HFAG values quoted above. 
The models with higher ml have a noticeably larger average value for u < 0.5 GeV, but 
the difference starts to become smaller for values of uj higher than this. This implies 
that raising ml tends to raise partial decay rates in the shape-function region, and that 
this effect is largest for the most restrictive cuts. We shall see in the next section that 
the changes in the shape-function model for different values of ml constitute the largest 
parametric uncertainty in the \Vub\ analysis. 

Since the emphasis of this paper is the study of the perturbative series for a given shape 
function model, in the numerical studies we limit ourselves to the shape-function model 
(17), with the moment constraints implemented as explained above. Recently, a different 
procedure for building shape-function models which satisfy the moment constraints was 
proposed in [27]. It would be interesting to see to what extent the numerical results 
change when using such models at common values of the HQET parameters and the scale 
yUo, but we do not explore this issue in the current work. 
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Figure 2: Dependence of the partial rate rl^"* (in units of |Vufepps~^) with a cut on 
P+ < A = 0.66 GeV on the matching scales /i^ and /Xj at LO (dotted), NLO (dashed), 
and NNLO (solid), with the parameter choices ml = 4.71 GeV and /i*^ = 0.2 GeV^. The 
case fih = f^i = fJ' corresponds to fixed-order perturbation theory (F.O.P.T.). 

3 Numerical results for leading-power partial rates 

In this section we study the effects of the NNLO perturbative corrections to the leading- 
power term For the time being, we limit ourselves to three benchmark partial rates, 
defined by the following cuts: 

1. P+ < 0.66 GeV 

2. Mx < 1.7 GeV 

3. El > 2.0 GeV 

Experimental data is available for each of these partial decay rates, and will be used 
to extract values of \Vub\ in the next section. The goal of this section is to address the 
following questions: 

• How does the NNLO analysis change the central values and error estimates compared 
to NLO? 

• Is resummation important, or is fixed-order perturbation theory sufficient? 
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Figure 3: Dependence of the partial rate rl^"* (in units of |Vu6pps^^) with a cut on 
Mx < Mq = 1.7 GeV on the matching scales fih and /ij at LO (dotted), NLO (dashed), 
and NNLO (solid), with the parameter choices ml = 4.71 GeV and /i*^ = 0.2 GeV^. The 
case fih = fJ'i = fJ' corresponds to fixed-order perturbation theory. 



• How well do the NNLO corrections in the large- /3o limit approximate the full results? 

We shall examine these issues in the following subsections, with the help of the information 
contained in Figures 2-7, and the numbers given in Section 3.3. Before interpreting the 
results, we first make some comments on the numerical evaluations. To study the behavior 
of the perturbative expansion, we give results where the structure functions (13) are 
evaluated at LO, NLO, and NNLO. In doing this, we evaluate the product of matching 
functions H ■ j and the RG-exponents S,r], a^j, a^i to the given order in a^. These 
quantities also depend on the renormalization scheme for m;,, which we choose as the 
shape-function scheme; the shifts from the pole scheme are made using the perturbative 
relations (16) as appropriate at a given order in a<j. At this point we evaluate the resulting 
expressions numerically, without any further expansion. This means, for example, that 
we do not use relations such as y'^r,Lo+6as ~ y°-r,ho[\ -|- fea^lny) (we have checked that 
the difference between the two expansions is negligible numerically). Finally, we use a 
common shape function when quoting results at LO, NLO, and NNLO, namely that where 
the moment constraints are implemented at NNLO. For reference, at ml = 4.71 GeV and 
= 0.2 GeV^, the shape-function model is S{iu) = 2.78 GeV'^n^-^^ exp {-il), where 
n = ujx 3.42 GeV-\ 
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Figure 4: Dependence of the partial rate Tu'^ (in units of |KifePps~^) with a cut on 
El > Eq = 2.0 GeV on the matching scales /i/i and /ij at LO (dotted), NLO (dashed), and 
NNLO (sohd), with the parameter choices ml = 4.71 GeV and /i*^ = 0.2 GeV^. The case 
fJ'h = fJ'i = fJ' corresponds to fixed-order perturbation theory. 



A further issue is the treatment of the charm-quark mass. The dependence on this 
parameter first enters at NNLO, through diagrams where a charm loop is inserted into a 
gluon propagator. In evaluating the perturbative expressions, we always set Uh = I and 
work with ri/ = 4 light flavors. In other words, we treat rric <^ fii and set rric = when the 
fermion loops needed for the calculation of the two-loop hard and jet functions contain 
a charm quark. The calculations from [23] show that this is not a bad approximation to 
the full dependence in the hard function H. Numerically, it may be more appropriate to 
treat rric ~ f^i, as discussed in the context of SCET in [28,29], but a consistent treatment 
would also require the rric dependence in the jet function, which is not yet available. 

A final complication is that the analytic expression for the hard function Hui contains 
a harmonic polylogarithm (HPL) of weight four, which cannot be expressed in terms of 
standard functions such as polylogarithms and their generalizations. For the numerical 
evaluation of HPLs we have used the Mathematica package HPL [30], and the FORTRAN 
subroutine hplog [31]. 
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Figure 5: Dependence of partial rates (in units of iKtbPps"^) on /ij, with fih = 
H'^/{1.5 GeV), at LO (dotted), NLO (dashed), and NNLO (sohd), for the cuts P+ < A = 
0.66 GeV, Mx < Mo = 1.7 GeV, and Ei > Eq = 2.0 GeV. In each case fx^^ = 0.2 GeV^. 



3.1 Impact of NNLO corrections on central values and errors 

The first question we wish to address is how the NNLO corrections and implementation 
of the dependence on the scale /Xj affect the central values and error estimates for the 
leading-power partial decay rates compared to the NLO analysis in [8]. We have shown 
in Figures 2-4 a set of plots meant to shed light on that issue. To compare central values 
obtained at different orders in perturbation theory we have shown in the upper-left hand 
plots of Figures 2-4 the dependence on the scale fih, for the fixed value fii = 1.5 GeV, 
which is the BLNP choice. We see that for this choice of fii the NNLO corrections are 
quite large in each case, and significantly lower the central values compared to NLO. 
At the same time, the perturbative uncertainties associated with the matching scale fih 
are reduced, as indicated by the flattening of the curves at higher orders in perturbation 
theory. 

A new element of our analysis compared to the procedure in [8] is that the matching 
scale fii is not necessarily fixed to the value /ij = 1.5 GeV. In the upper right-hand plots 
of Figures 2-4, we display the dependence of the results on the scale /ij, for fixed values 
Hh = iTT'l/'^- note a large dependence on the intermediate scale at NLO, which 
is reduced but still significant at NNLO. Because the dependence on fii is so strong at 
NLO, even small changes of the value of fii can can alter rather drastically the agreement 
between the NLO and NNLO results. To illustrate this effect we have shown in the bottom 
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Figure 6: Dependence of partial rates (in units of |V^fepps~^) on ml at LO (dotted), NLO 
(dashed), and NNLO (solid), for the cuts P+ < A = 0.66 GeV, Mx < Mq = 1.7 GeV, 
and Ei> Eo = 2.0 GeV. In each case fxf = 0.2 GeV^ 



left-hand plots of Figures 2-4 the dependence of the partial rates on /x/j for /Xj = 2.5 GeV, 
which is higher than the BLNP choice. This higher value of the intermediate scale brings 
the NLO results into closer agreement with the NNLO results. 

So far, we have fixed one of the two matching scales to a default value and studied 
the behavior of the partial decay rates under variations of the other. Another option is 
to establish a correlation between the scales and vary them simultaneously. Setting them 
equal to one another corresponds to fixed-order perturbation theory and will be discussed 
in the next section. A different choice, which is somewhat more natural in SCET, is to 
vary the scales such that the relation fth ~ is respected, where fig is a typical soft 

scale. In Figure 5, we show the behavior of the partial rates under such a correlated 
variation for the choice fig = 1-5 GeV. It is remarkable that the NLO result is almost 
entirely stable under scale variations. This feature appears to be accidental, since at 
NNLO the scale dependence actually becomes stronger under this correlated running. 

In addition to the unphysical dependence on the matching scales fih and /ij, the partial 
rates contain a rather strong parametric dependence on the numerical value of the b- 
quark mass ml, mainly through the moment constraints on the shape-function model; 
scaling relations for this dependence were derived in [8]. We have already shown in 
Figure 1 how changing ml distorts the shape-function model. We show in Figure 6 
how these changes in the model translate into changes in the partial rates in the range 
4.65 GeV < ml < 4.77 GeV. Obviously, the present uncertainty in the value of the 6-quark 
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Figure 7: Comparison of the NNLO partial rates rl°^ (in units of |V^u6pps~^) in resummed 
perturbation theory and fixed order, for the cuts P+ < A = 0.66 GeV, Mx < Mq = 
1.7 GeV, and Ei > Eq = 2.0 GeV. The dashed hne is the variation of the resummed 
results with /i/j = (4.25 GeV)x, for /ij = 2 GeV, the dotted line is the variation with 
/ij = (2 GeV)x, for fi^ = 4.25 GeV, and the solid line is the variation of the fixed order 
result with ^={3 GeV)x\ In all cases m* = 4.71 GeV and fif = 0.2 GeV^. 



mass adds a large parametric uncertainty to the analysis. 

From the above discussion we can conclude that, compared to BLNP analysis in [8], 
where /ij = 1.5 GeV, the net effect of the NNLO corrections is to lower the central values 
for partial rates by around 15-20%, while at the same time reducing the perturbative 
uncertainty on the scale fit- Moreover, we have pointed out that there is a considerable 
dependence on the scale /ij at NLO, which is reduced but still significant even at NNLO. 
To make these statements more quantitative, we give central values and error estimates 
for some input parameters. In particular, we choose the values iJ,f'^^ = 2.0 GeV and 
f^'h^ = 4.25 GeV, and then vary the scales in the range l/\/2 < ^ih^i/p^'^i < V^- We then 
find for the three benchmark partial rates, in units of |K(6pps~^: 

• P+< 0.66 GeV: 





J- u 






NLO 
NNLO 


60.37 
52.92 


+3.52 
-3.37 

+ 1.46 
-1.72 


+3.81 
-6.67 

+0.09 
-2.79 
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• My < 1.7 GeV: 





J- u 






NLO 
NNLO 


61.86 
52.26 


+3.23 
-3.21 

+ 1.05 
-1.40 


+ 1.89 
-5.50 

+0.65 
-3.90 


El > 2.0 GeV: 




r(o) 
J- « 


jJ'h 




NLO 
NNLO 


25.98 
22.61 


+ 1.63 
-1.61 

+0.63 
-0.78 


+ 1.69 
-2.81 

+0.01 

-0.94 



The default numbers refer to the value of the partial rate for /i/^ j = /x^*^/. To obtain the 
uncertainties associated with variations of fih, we keep yUj fixed to its default value, and 
then assign upper (lower) errors by picking out the highest (lowest) value of the decay 
rate in the range 1/V2 < jjt!^^ < the procedure for /x, is analogous. 

3.2 Comparison with fixed-order perturbation theory 

In this section we compare the results from RG-improved perturbation theory with those 
obtained in a fixed-order calculation. The standard argument to motivate the use of 
resummed perturbation theory is that even though logs of are not particularly 

large, the running coupling constant runs quickly at low scales /ij ~ 1.5 GeV, so the 
scale separation between the hard and jet scales is important. If this were the case, we 
would expect to see a large scale dependence in the fixed-order result. To test whether 
this actually happens, we have shown in the bottom right plots of Figures 2-4 the results 
where the matching scales are set to fih = f^i = Z^, which corresponds to fixed-order 
perturbation theory. We see that the NNLO results are rather stable under variations of 
fi from 1.5 GeV to 2m^. Also, the convergence of the perturbation series is actually better 
than in resummed perturbation theory, in the sense that the NNLO results lie within the 
ranges covered by varying the renormalization scale /x in the LO and NLO results. 

To illustrate further the differences between resummed and fixed-order perturbation 
theory we show in Figure 7 the NNLO results for both cases. After the substitutions 
indicated in the caption, the curves can be used to compare the fixed-order results in the 
range /i = (1.5 — 6) GeV with the resummed results in the ranges 1/ V2 < fih,i/ < 
with /if'^^ = 2.0 GeV and = 4.25 GeV. In this way, the scale in the fixed-order result 
is varied in a range which very nearly covers the lowest value of /Xj to the highest value 
of fih- To compare results more closely, we list tables of central values and errors, to be 
compared with the resummed results in the previous section. We choose the default scale 
as /i = 3.0 GeV, and quote the uncertainty obtained by varying it up and down by a 
factor of two. We then find, in units of |V„fepps~^: 

• P+< 0.66 GeV: 



Fixed-Order 


p(0) 
J- u 


^^ 


NLO 
NNLO 


49.11 
49.53 


+5.43 
-9.41 

+0.13 
-4.01 
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• My < 1.7 GeV: 



Fixed-Order 


J- u 


/i 


NLO 


51.81 


+3.69 
-8.62 


NNLO 


50.47 


+0.01 
-2.62 


El > 2.0 GeV: 


Fixed-Order 


r(o) 

J- u 


/i 


NLO 


21.01 


+2.04 

-3.54 


NNLO 


20.99 


+0.04 
-1.43 



At NNLO the errors are comparable with the resummed results from the previous section. 
In each case the central values for the fixed-order results are noticeably lower than the 
resummed ones. 

The analysis above indicates that, for the leading-order term in the SCET expansion, 
the factorization of the perturbative coefficient multiplying the leading-order shape func- 
tion into jet and hard functions is not strictly necessary: using the fixed-order results does 
not lead to large scale uncertainties compared to the treatment in RG-improved pertur- 
bation theory, nor to a poor convergence of the perturbative expansion. Given this fact, 
it is legitimate to ask whether "kinematic" power corrections, namely those suppressed 
by the perturbative ratio of jet to hard scales and scaling as P+/P_ ~ Aqcd/ rrib in the 
shape-function region, should be treated separately from the leading-order term. The one- 
loop analysis in [8] shows that the 1 / expansion of these kinematic power corrections 
converges very quickly, and that the leading-order term is indeed of the size expected of a 
l/rrLh correction. To check this at two loops would require the calculation of the hadronic 
tensor in fixed-order perturbation theory to this order, which has yet to be done. While 
partial results in the large-/3o limit are available, we shall see in the next section that, for 
the leading-order term, they do not accurately approximate the full two-loop calculations. 
In absence of evidence to the contrary, we shall assume that ignoring the kinematically 
suppressed two-loop corrections provides an accurate approximation of decay rates to this 
order. 

3.3 NNLO corrections in the large-/3o limit 

We now give illustrative results for the case where the large-/3o approximation to the 
two- loop hard and jet functions is used. To obtain these results, we pick out the niCp 
color structure in the corrections to H ■ j, and then set ri; — 3/3o/2. As a check on 
our results, we have confirmed that convoluting the resulting expression with the large-/9o 
limit of the partonic shape function evaluated at NNLO in [26] , and setting all matching 
scales to a common scale we reproduce the NNLO corrections to the partial rate with 
a cut on P+ obtained in the large-/3o limit in [32] (these can also be obtained from the 
triple differential result in [33]). To compare numerical values of partial decay rates in the 
large-/3o limit with the exact ones, we evaluate the hard and jet functions at a common 
scale /i/i = yUj = 1.5 GeV. Furthermore, we choose ml = 4.71 GeV and yU*^ = 0.2 GeV^, 
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and convolute both results with a common shape function obtained by implementing the 
moment constraints at NNLO. Then the full NNLO results for labeled "QCD", and 
the large-/3o results, labeled "BLM", read for the three cuts (in units of \V^%\ ps~^): 



P+ < 0.66 GeV: 

QCD: 54.43 + 0.11 [a, 
BLM: 54.43 + 0.11 [a, 

Mx <1.7 GeV: 

QCD: 58.06 - 2.76 [a, 
BLM: 58.06 - 2.76 [a, 

El > 2.0 GeV: 

QCD: 22.85 + 0.21 [a, 
BLM: 22.85 + 0.21 [a, 



+ ( - 3.68 [h] - 0.26 [j] - 4.61 [hj] = -8.55) 
+ ( - 14.1 [h] + 14.1 [j] = -0.02) 

+ ( - 4.05 [h] - 1.99 [j] - 4.04 [hj] = -lO.l) 
+ ( - 15.7 [h] + 13.0 [j] = -2.68) 

+ ( - 1.64 [h] + 0.65 [j] - 2.28 [hj] = -3.27) 
+ (- 7.12 [/i] +6.32 [j] =-0.80) 



a„ 



45.99 
54.52 

45.22 
52.62 

19.79 
22.26 



We have included labels to distinguish the contributions from different orders in a^, and 
for the NNLO pieces we have also indicated whether the contribution comes from the hard 
function {[h]), the jet function {[j]), or the product of one-loop hard and jet functions 
([^j]). We see that the contributions of the two- loop hard and jet functions in the 
large-/3o approximation are larger than the full results. In each case, they undergo large 
cancellations in the sum and give results which are much smaller in magnitude than the 
complete ones. Given the poor agreement of results in the large-/9o approximation with 
the full NNLO results for the leading order term in the HQET expansion, it is somewhat 
questionable that including only such terms leads to any phenomenological improvement 
compared to NLO. For the jet function j, the poor agreement between the full results and 
those in the large-/3o limit was previously noted in [24]. 



4 Impact on the extraction of \Vub\ 

In this section we study the numerical impact of our results on the extraction of [Vub\ 
from experimental information on partial decay rates. In order to do so, we combine our 
NNLO results for the leading-power term ri°^ with the power corrections obtained in [8]. 
A given partial rate is then obtained in the form 

Tn = fW + [(Fl^-(i) + Tl^''^'^) + (Fl^^-^^^) + Tl^''^'^)] . (18) 

The power corrections from the BLNP analysis are split into hadronic contributions, 
involving subleading shape functions [34-36] , and kinematic contributions, which account 
for terms suppressed by powers of P+/Mb ~ P+/P- and thus scaling as Aqcd/-^s in the 
shape-function region^. The hadronic terms are treated at tree level, and the kinematic 
terms at one loop. 

^For a recent analysis of these corrections in terms of subleading jet functions in SCET. see [37]. 
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We summarize our results in the tables of numbers below, for the three partial rates 
studied in the previous section, as well as for different values of the Ei and Mx cuts. 
The central values and errors are obtained as follows. First, we evaluate the leading- 
power term Tu^ as in Section 3.1. We choose the default scales as nf^^^ = 4.25 GeV and 
/if ^ = 2.0 GeV, and the HQET parameters as m* = 4.707 GeV and fx*^^ = 0.216 GeV^. 
Then, we evaluate the power-suppressed terms exactly as in [8], for the same choice of 
HQET parameters as for the leading term. Adding these two numbers together gives a 
default value for F^. To this default value we assign uncertainties coming from a number 
of different sources. To associate uncertainties with the HQET parameters, we vary them 
simultaneously in the leading-order term and power corrections, and assign errors for the 
ranges m* = (4.707lEJ:[J|) GeV and fif = (0.216l[5:|]^^) GeV^. To estimate perturbative 
uncertainties associated with the choice of matching scales in the leading term, we vary 
the default choices for jj^ and Hi up and down by a factor of \/2, and add these in 
quadrature to obtain the uncertainty labeled fih,i in the tables. Perturbative scales also 
appear in the power corrections, through resummation factors, and through the scale /2 
for the kinematic power corrections. We vary these as in [8] , and add the different sources 
of perturbative uncertainty in quadrature to obtain the total perturbative uncertainty for 
the power corrections, labeled /Xpow in the tables. Finally, we associate errors with the 
functional variation of the subleading shape functions, called SSF in the tables, and add 
an uncertainty of ST^-^ = ±1.26\Vub\'^ps~^ to take account weak annihilation effects (not 
shown in the table). The analysis shows that the uncertainty from ml is the dominant 
one. To quote the uncertainties in such a way that isolates this effect, we add all of 
the others in quadrature to obtain a total uncertainty excluding that associated with ml, 
which we instead quote in the last column of the tables. The results read (in units of 
iKfePps-^): 

• P+ < 0.66 GeV: 







l^h,i 


/^pow 




SSF 


tot 


ml 


NLO 
NNLO 


54.28 
46.77 


+5.07 
-7.36 

+ 1.44 
-3.29 


+ 1.91 
-1.48 

+ 1.91 
-1.48 


+ 1.27 
-1.07 

+2.08 
-1.56 


±1.41 
±1.41 


+5.87 
-7.81 

+3.69 
-4.36 


+8.20 
-7.05 

+6.34 
-5.48 


Mx < 1.55 GeV: 




F„ 


f^h,i 




f^f 


SSF 


tot 


ml 


NLO 
NNLO 


54.73 
47.09 


+4.68 
-7.08 

+ 1.22 
-3.44 


+ 1.99 
-1.31 

+ 1.99 
-1.31 


+ 1.02 
-0.94 

+ 1.79 
-1.42 


±1.30 
±1.30 


+5.50 
-7.48 

+3.46 
-4.34 


+8.02 
-6.87 

+6.26 
-5.41 


Mx < 1.7 GeV: 




r„ 


l^h,i 


/^pow 


l^f 


SSF 


tot 


* 


NLO 
NNLO 


61.21 
51.60 


+3.69 
-6.29 

+1.24 
-4.13 


+2.48 
-1.60 

+2.48 
-1.60 


+ 1.85 
-1.46 

+2.28 
-1.72 


±0.63 
±0.63 


+5.02 
-6.80 

+3.85 
-4.95 


+ 7.36 
-6.42 

+5.80 
-5.05 
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Method 


AS«-p [10-4] 


Kb\ [10-^] 

NLO 




Kbl [10-=^] 

NNLO 


El > 2.1 GeV 


3.3 ±0.2 ±0.7 


3 56 ± AQ+O-^s+O-Si 

O.OU IE U.^U_Q 27_o.26 


3.81±0.43+[Jf+°iJ 


CLEO [38] 










El > 2.0 GeV 


5.7 ±0.4 ±0.5 


0.9 { ± O.ZZIq 


H0.26 
-0.25 


4 30 ± 24+°-26+0-28 


BABAR [39] 










El > 1.9 GeV 


8.5 ± 0.4 ± 1.5 


4.27 ±0.39+°;?^: 


1-0.25 
-0.22 


4.65 ± 0.43l°:L^l°JI 


BELLE [40] 










Mx < 1.7 GeV 


12.3 ± 1.1 ± 1.2 


3.55±0.24+°;22: 


1-0.21 
-0.19 


3.87 ±0.261™:- 


BELLE [41] 










Mx < 1.55 GeV 


11.7 ±0.9 ±0.7 


3.67±0.18lE];??: 


hO.26 
-0.24 


3.96±0.19l°i°l°:^^ 


BABAR [42] 










P+ < 0.66 GeV 


11.0 ±1.0 ±1.6 


3.56±0.3i;°:??^°:i 


3.84 ± 0.331°:?^!°:^^ 


BELLE [41] 










P+ < 0.66 GeV 


9.4 ± 1.0 ±0.8 


3.30 ± o.23+!J:^^: 


f-0.25 
-0.22 


3.55 ± oMtritr.t 


BABAR [42] 











Table 1: Values of \Vub\ determined at NLO and NNLO, for the parameter values discussed 
in the text. The uncertainties in the experimental measurements of AS'^^p are statistical 
and systematic, respectively. In the columns labeled \Vub\ the first error is experimental, 
the second is the sum of all theoretical and parametric errors except for that from ml, 
and the third is that from m^. We have combined errors of the same type by adding in 
quadrature, and used tb = 1.584 ps for the average 5-meson lifetime. 



• Ei> 1.9 GeV: 









/^pow 




SSF 


tot 


ml 


NLO 
NNLO 


29.21 
24.65 


+ 2.14: 

-3.29 

+0.63 
-1.61 


+ 1.52 
-1.01 

+ 1.52 
-1.01 


+0.63 
-0.51 

+0.75 
-0.60 


±0.53 
±0.53 


+3.03 
-3.73 

+2.27 
-2.41 


+3.51 
-2.91 

+2.95 
-2.43 


El > 2.0 GeV: 




r„ 


f^h,i 




f^f 


SSF 


tot 


* 


NLO 
NNLO 


22.76 
19.38 


+2.34 
-3.23 

+0.63 
-1.22 


+ 1.19 
-0.88 

+ 1.19 
-0.88 


+0.44 
-0.35 

+0.60 
-0.46 


±0.59 
±0.59 


+3.00 
-3.64 

+2.03 
-2.10 


+3.23 
-2.66 

+2.70 
-2.21 


El > 2.1 GeV: 




r„ 


l^h,i 


f^pow 




SSF 


tot 


ml 


NLO 
NNLO 


16.30 
14.17 


+2.50 
-3.05 

+0.64 
-1.06 


+0.97 
-0.89 

+0.97 
-0.89 


+0.19 
-0.16 

+0.38 
-0.29 


±0.70 
±0.70 


+3.05 
-3.49 

+ 1.89 
-2.02 


+2.87 
-2.31 

+2.39 
-1.93 
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From these numbers and the corresponding experimental results, we arrive at the values 
of \Vub\ listed in Table 1. An examination of the table shows that the NNLO corrections 
shift the values of \Vub\ upwards by roughly 10% compared to NLO. 

5 Conclusions 

We studied the impact of NNLO perturbative corrections on the leading term in the 
l/rrih expansion for partial decay rates in B ^ X^Wi decays. These corrections were 
implemented within a modified form of the BLNP framework, which allows for variations 
of both the perturbative jet scale and the hard-matching scale fih in the resummed 
partial rates. The particular choice /ij = fit corresponds to fixed-order perturbation 
theory, which allowed us to perform a detailed comparison between fixed-order results and 
the resummed results from BLNP. Within resummed perturbation theory, we found that 
the dependence on the intermediate scale fii introduces sizeable perturbative uncertainties 
at NLO, which are reduced but still significant at NNLO. For the conventional choice 
Hi = 1.5 GeV used in previous analyses in the BLNP formalism, the NNLO corrections 
also induce fairly large downward shifts in the partial rates. For higher values of and 
in fixed-order perturbation theory, these shifts are more moderate. We also compared 
between the full NNLO results in fixed-order perturbation theory and those obtained 
in the large-/3o approximation. We found that the large-/5o approximation for the NNLO 
corrections provides a poor approximation to the full results, at least for the leading-order 
term in the 1 / rrih expansion. Whether this would also be true upon the inclusion of terms 
suppressed by kinematic factors of P+/Mb-, which are small in the shape-function region, 
is an open question. Finally, we combined our new results for the leading-order partial 
rates with the known power corrections up to 1/m^, and showed how our analysis impacts 
the determination of \Vub\ from several experimental measurements. For parameter and 
scale choices typically used in current analyses within the BLNP framework, the effect of 
the NNLO corrections is to raise the value of \Vub\ by slightly less than 10% compared to 
NLO; the exact results are shown in Table 1. 

Acknowledgments: C.G. is partially supported by the Swiss National Foundation as 
well as EC-Contract MRTN-CT-2006-035482 (FLAVIAnet). The Albert Einstein Center 
for Fundamental Physics (Bern) is supported by the "Innovations- und Kooperationspro- 
jekt C-13 of the Schweizerische Universitatskonferenz SUK/CRUS". 
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A Appendix 

A.l The hard and jet functions 

The hard functions Hui are derived from the SCET Wilson coefficients Cj, which arise 
when matching the QCD b ^ u transition current onto SCET. The two-loop QCD cal- 
culations needed for the NNLO analysis were recently completed in [20-23]. In terms of 
the quantities = CiCj, the Hui read 

Hui = Hii, Hu2 = 0, Hu3 = h Hi2 + H23 + 7-f^22 • (Al) 

y 4 

The Laplace transformed jet function j was calculated to NNLO in [24]. To two-loop 
order, the result is 



An 



-3L + 7 



o - 
ZTT" 



a \ 2 

-Cf[^) [CFJF + CAJA + TpnfJf], (A2) 



where 



Jf = 2L* -61-" + { — 



37 47r2 



45 



9 V 18 3 y 

53129 1557r2 STn^ 



y +47r2-24C3)i. + ^ 

3155 llvr^ , , . , 
+ ^ + 40C3U 



20 5 9 77r2 Glvr^ 



12 



54 



9 



648 



Jf 



9 9 



36 180 
494 47r2 



18C; 



3 ) 



L 



4057 137r2 



162 



9 



90 



-6C3 



(A3) 



A. 2 RenormaUzation-group factors and anomalous dimensions 

Here we list the perturbative expansion of the renormalization-group functions in (10) 
up to NNLO. To do this, we first define the expansion coefficients of the cusp anomalous 
dimension and QCD /3-function as 



rcusp('-^s) 



in/ 



An J 



-2a. 



An) 



a. 



(A4) 



and similarly for the other anomalous dimensions. In terms of these quantities, the 
function ar is given by [17, 19] 



+ 
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Ji 
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To 




Po 




Po 



327r2 



+ . . . 



(A5) 
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Note that this resuit invoives the three-ioop anomaious dimension, which is known for the 
cusp anomalous dimension and for a^j, but not for the anomalous dimension a^i. The 
result for the Sudakov factor S is [19] 



An 



Inr 

r 



r + In r) + In^ r 



+ 



Air 



[1 — r + rlnr) + 



PI /3o 



[1 — r) Inr 



/3o /3o^o To 



'1-r) 



+ 



Ati 



2/3o /3o Uo 



£3 _/53 

ro /9o 

3^_£3 

4/3o 



2^ , /5? 



/5o 



+ 



2r, /?,r_\ ]^ 



^ ATiX r2 
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P2^X 



Inr 



2ro A' 
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:i-r) 



r 



(A6) 



where r = as(/x)/Q!s(i^). Whereas the three-loop anomalous dimensions and /3-function 
are required in (A5), the expression for S also involves the four- loop coefficients Fa and 
Ps- 

We now list expressions for the anomalous dimensions and the QCD /^-function, quot- 
ing all results in the MS renormalization scheme. The expansion of the cusp anomalous 
dimension Fcusp to two- loop order was obtained some time ago [43], while recently the 
three-loop coefficient has been obtained in [44]. For the four-loop coefficient F3, we fol- 
low [45] and use its [1,1] Fade approximant, F3 ^ Fg/Fi. The results are 



To = ^Cf 
Fi = ACf 

To = ACf 



245 



1347r2 llTT^ 22 \ ^ ^ / 418 AOn^ 



56 



+ CfTfUj [ -y 



16Cs 



16^2 2 

27 ^ ^ 



7849, 4313, 1553 for Uf = 3, 4, 5 



(AT) 



The SCET anomalous dimension 7' was deduced at two loops using RG-invariance 
along with results for the jet and shape function anomalous dimensions in [17], and was 
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confirmed through the exphcit calculations in [20-23] . The result is 

7o = -5C^, (A8) 



„^ r/3 ^ \ ^ /1549 Ttt^ 11 , \ ^ /125 tt^\ 

7 = -SCf + 3C3 + + Cs Ca- { + — ]nf 

[\16 4 ^7 V 432 48 4 ^7 V216 24/ ^ 

To evaluate the RG-factor ay at NNLO requires the three-loop anomalous dimension, 
which is not yet known. For that, we use the [1,1] Fade approximation, 72 ~ 7f/7o- 
Although this same approximation works well for F^usp it works poorly for 7j, as can be 
verified using the explicit results listed in [19]: 



1769 Il7r2 ^ ^ /242 4tt^ 





— — 3C p 








li 


= cl[- 


- + 27r2 - 
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24C3) 


+ CfCa 
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Stt^ 
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- 68C3 + 



54 9 '7 ^ V 27 9 
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C3 + 24OC5 

+ C|C.(-i^ + ^ + ^-^C3-^C3-120C. 
F A\ ^ 9 135 3 ^3^ ^ 

2 / 412907 4197r2 197r^ 5500 SStt^ 

+ O pO 4 (,3 (,•? — Z0ZC5 

^ ^ 2916 243 10 9 ^9 ^ ^ 

+ ClTpUf \ C3 

^ ^ V 27 9 135 9 ^ 

^ ^ ^ / 5476 llSOvr^ 467r^ 2656 , 

^ ^2 2 / 13828 807r2 256 \ , 

+ CpTWf C3 • A9 

^ ^ V 729 81 27 ^7 ^ ^ 

Finally, the expansion coefficients for the QCD /3-function to four-loop order are 

^0 = y C'a - ^TpUf, 
34 20 

A = y - y CaTfUj - ACfTfUj , (AlO) 

ft . 2|I ci . - f - IIH ci) r.„, + (f + f c.) t|„5 

^■^6 ^ V 162 27 ^7 ^ V 162 81 ^7 ^ 729 ^ 

where jS^ is taken from [46] and is evaluated for Nc = 3. 
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